Dualities for Multi-State Probabilistic Cellular 
Automata 

o 

F. Javier Lopez 1 , Gerardo Sanz 1 and Marcelo Sobottka 2 

1 Departamento de Metodos Estadfsticos and BIFI, Facultad de Ciencias, Universidad 
de Zaragoza, C/ Pedro Cerbuna 12, 50009 Zaragoza, Spain. 
, 2 Dcpartamento de Ingenieria Matematica, Facultad de Ciencias Fisicas y 

Matematicas, Universidad de Concepcion, Casilla 160-C, Correo 3, Concepcion, Chile. 



CM 

c3 



E-mail: j avier . lopez@unizar . es, gerardo . sanzQunizar . es, msobottkaOudec . cl 



Abstract. 

In this paper a new form of duality for probabilistic cellular automata (PCA) is 
introduced. Using this duality, an crgodicity result for processes having a dual is 
proved. Also, conditions on the probabilities defining the evolution of the processes for 
the existence of a dual arc provided. The results are applied to wide classes of PCA 
which include multi-opinion voter models, competition models and the Domany-Kinzcl 
, ■ model. 

o 
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^ 1 1. Introduction 

^— > 

a 

Probabilistic Cellular Automata (PCA) are discrete-time stochastic processes with state 
j> . space X := W zd , where W = {1, . . . , M} and d > 1, with finite range interactions on 
^ , the integer lattice Z d (this is interpreted as at each point z G Z d there is a particle which 
take values in the set W; the points x G X are called configurations). More precisely, 
suppose N C Z, d is a finite subset and / : W N x W — * [0, 1] is a transition function, 
then a PCA is a discrete-time homogenous Markov process rj t = {f]t(z) G W : z G Z d } 
whose evolution is given, for all s G N, z G Z d , w G W, and x G X, by 

P {t] s +i{z) =w\r] s = x} = f{{x{z + n)) neN , w). (1) 

In other words, t] t is an interacting particle system on X in discrete time such that the 
transitions at different sites are conditionally independent, given the current state of the 
system. For definitions and main results on PCA, the reader can consult [1] and [2]. 

The set X = W zd is endowed with the product topology. Let "P(X) be the set of 
probability measures on X with the weak* convergence topology. Note that "P(X) is 
compact with respect to this topology. 

Let rjt and £t be two stochastic processes (which can evolve either in continuous or 
discrete time) with state spaces X and Y respectively, and suppose H : X x Y — ► R is 
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a bounded measurable function, then r\ t and £ t are dual to one another with respect to 
H if for all x G X and y G Y (Definition II.3.1 of [3]) 

E V0=x [H(r) S: y)]=-Es 0=y [H(x : Q] : (2) 

where E 70=7 [-] denotes the expectation when the process j t starts with configuration 7. 

If rjt and £ t are both discrete-time Markov chains with transition matrices P and 
Q respectively, then the previous expression for duality can be written as 

P S H = H{Q S ) T , (3) 

where the superscript T stands for the transpose. Furthermore, due to the Markovian 
property, if the duality equation holds for s = 1 then it holds for all s G N. 

Duality allows obtaining relevant information about the evolution of r] t , which 
usually has uncountable state space, by studying the evolution of £t which in general is 
chosen having countable state space. In fact, duality has been widely used in the study 
of spin systems (continuous-time interacting particle systems with W = {0, 1}), where 
Y is the set of all finite subsets of Z d (see [3], [4], [5], and [6]). For the case where rj t is 
a continuous-time interacting particle system with \W\ > 2, Lopez and Sanz [7] studied 
an extended version of (2), which is given by the equation: 

E V0=X [H(r) s ,y)} = E^ 0=y [h (x, f s )e"'o v ^ du j , (4) 

for all x G X and y G Y, and for some function V : Y — > [0, 00). In particular, this 
duality equation was used in [7] to obtain results about the long time behaviour of 
the continuous-time multi-opinion noisy voter model and the continuous-time 3-opinion 
noisy biased voter model. 



For the discrete-time case, Katori et al. [8] develop a theory of duality for two state 
PCA based on the equation 

c l^nB|l j (5) 

where c is a parameter depending on the processes t] t and £ t ; A,B C Z, with at least 
one of them being finite; 77^ and £,f represent the set of particles with value 1 at time s 
starting from configurations where only the particles in A and B, respectively, have the 
value 1; and | • | stands for the cardinality of a set (this equation was first considered 
for interacting particle systems by Sudbury and Lloyd in [9]). In [8], this relationship 
is used to study duality for some cases of the Domany-Kinzel model, which is the PCA 
with state space {0, 1} Z , and evolution defined by the parameters ao, 01,02 G [0,1] as 
follows: 

F{rj s+1 (z) = 1| r] s (z - 1) = w_ u r) s {z) = w ,rj s (z + 1) = w t } = a w _ 1+Wl (6) 

When ao = and a<i < 01, Katori et al. [8] have presented a dual for r\ t . Such 
dual is a model in which each stage is updated by a ai-thinning of all sites, followed 
by an application of the Domany-Kinzel model with parameters a' Q — 0, a' x — 1, and 
a' 2 = a%ja\. When a\ > 02, they also studied the limit behaviour of the process. The 
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technique used in [8] was to locate the process rj t in a finite subset of Z to represent it 
by a transition matrix P, and so to find H and Q which satisfy the duality equation (3). 
The technique used by Katori et al. depends strongly on the fact that ao = (see also 
chapter 5 of [4] for results on this model, with a = 0, using percolation theory). In [10] 
and [11], those results are extended to a larger class of one dimensional two state PCA. 

For multi-state PCA, an equation analogous to (5) is introduced in [12] to define 
self-dual processes (that is, processes r\ t that satisfy the equation with £ t = rj t ) and give 
conditions for self-duality 

Our objective in this paper is to give a general theory of duality for multi-state PCA. 
The duality equation we define (see (7) below) is given for general duality functions H 
and is not restricted to self-duality. Moreover, the introduction of an auxiliary function d 
in the equation allows the inclusion of wider classes of processes for each dual function H. 
We will show the usefulness of our concept by obtaining an ergodicity result (Theorem 
2.1) for processes satisfying (7). Equation (7) and Theorem 2.1 are stated for general 
H and we show their applicability with two particular cases of H which cover two 
wide classes of PCA, providing conditions on their transition probabilities such that 
they satisfy (7) and, in that case, providing conditions for their ergodicity (that is, the 
existence of an unique equilibrium measure and the convergence of the process to that 
measure). We will also apply our results to some examples. 

The duality equation we propose is 

Definition 1.1. Given two discrete-time Markov processes, rjt with state space X and 
£t with state space Y, and H : X x Y — > R and d : Y — *■ [0, oo) bounded measurable 
functions, we say rjt and £ 4 are dual to one another with respect to (H, d) if 

E vo=x [H( Vl , y)\ = d(y)E^ =y [H(x, &)] . (7) 

Notice that the duality equation with respect to (if, d) is only written for the one- 
step evolution of the processes rj t and £ t as they are both Markov processes, but an 
induction argument yields the duality equation for the s step evolution of the process: 

E V0=x [H( Vs ,y)} = d(y)E $0=v [d^ 1 ) ■ • • 6)]- (8) 

When rjt and £ t are Markov chains with transition matrix P and Q, respectively, 
equation (7) can be written as 

PH = H(DQ) T , (9) 

where D is the diagonal matrix with D yy = d(y). In this case, the general expression 
for time s > 1 is P S H = H ((DQ) S ) T . 

The advantage of considering the function d on the right side of (7) is that processes 
that do not have dual with respect to some function H in (2) may have it with respect 
to (H, d) (see Remarks 3.9 and 3.15 below). When d = 1, we have the classical notion of 
duality (2), which means that the evolution of r\ t is easily understood from the evolution 
of £t. When d ^ 1, the relationship between rjt and £t becomes more complicated; 
nevertheless, in that case, important information about rjt can be obtained from £ t . 
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Although equations (7) and (8) can be seen as a discrete counterpart of equation (4), 
given for the continuous-time case, there are several differences between the techniques 
used to exploit the former and the latter. In fact, instead of using directly equation (8) 
to prove results on ergodicity for the process, we show that there is a transformation £ t 
of the dual process such that r] t and £t are dual in the classical sense (with respect to a 
modified function H), see Theorem 2.1(a) below. This approach is useful, for instance, 
for writing down the explicit form of the invariant measure (11) and for providing 
alternative conditions for ergodicity (see Remark 2.2) which are used, for instance, in 
the study of the Domany-Kinzel model (Corollary 3.8). 

The organization of the paper is as follows. In §2 we present the basic tools for the 
duality of stochastic processes and obtain conditions on H and d for the ergodicity of 
PCA. In §3 we analyze, using the duality theory developed, two general classes of PCA. 
For each of these classes, we first obtain conditions on the transition probabilities of 
the processes so that they have a dual and then give conditions for their egodicity. The 
results are applied to specific models such as the multi-opinion noisy voter model, the 
Domany-Kinzel model and competition models. For the particular case of the Domany- 
Kinzel model with parameters < oq < a\ < 02 we find a dual process for many cases 
not covered by [4] or [8]; and we get new results for the ergodicity of this model. 

2. Dualities 

In this section we consider probabilistic cellular automata rjt for which there exist a pair 
of functions (H,d) and a dual process £t with respect to (H,d). Sufficient conditions 
on (H, d) for the ergodicity shall be presented together with the characterization of the 
corresponding equilibrium measure. 

In what follows, we first recall some basic definitions and results about dualities for 
stochastic processes. Suppose r\ t on X and £ 4 on Y are two stochastic processes that are 
dual to one another with respect to a function H . Given a probability measure [i on X, 
define for any i/GY 



For any s > 0, denote by jj, s the distribution obtained from the initial distribution \x 
when r] t evolves until time s. In our setting, X = W zd and rjt is a PCA given by (1); then 
fi s is defined on the cylinder w = [w Zl , . . . , w Zm ] := {x G X : x(zi) = w Zi , 1 < i < m} 



uec s 

where C s is the family of all cylinders of X defined on the coordinates 
{zi + rij : i = 1, . . . , m; rij G N}. 




by 
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Therefore, we have that 



Mv) = / H(x,y)dfi s (x) = / E m=x [H(rj s ,y)]dn(x 



and by the duality equation (2) it follows that 

Mv) = I^ E v =x[H(r] s ,y)]dfi{x) 

= J x ^o=y[ H ( x ,ts)}d[i(x) 

(10) 

= %=y lIx H ( x ^«) d t M ( x )] 

= %=« [£(6)] ■ 

Suppose the set of the linear combinations of the functions {H(-,y) : y G Y} is 
dense in C(X), the space of continuous real functions on X. If for any y, E^ 0=y [fi(C, s )} 
converges to some function of y (not depending on /x) as s goes to oo, then for any 
/ G C(X), it follows that J* /d// s converges to some functional £(/), which means /i s 
converges in the weak* topology to some probability measure v. Furthermore, if the 
above holds for any initial distribution /i, it implies the ergodicity of r\ t : that is, the 
convergence in distribution of the process to the equilibrium measure v from any initial 
measure. 

Now we are able to study the case of the duality with respect to (H, d). 



Theorem 2.1. Suppose r\ t is a Markov process with state space X and £ t is a Markov 
chain with state space Y, which are dual to one another with respect to (H,d), where 
< d(y) < 1 for ally G Y. Then: 

(a) There exist a Markov chain £ t with state space Y = Y U {p}, where p is a new 
state, and a bounded measurable function H : X x Y — > R such that rj t and £t are 
dual to one another with respect to H; 

(b) denoting by the set of all absorbing states of £t, if 

(i) the linear combinations of {H(-,y) : y G Y} is a dense set o/C(X); 

(ii) d(y) < 1 for any y ^ ; and sup d(y) < 1; 

y&Y: d(y)<l 

(in) H(-, 9) = c{6) for allO G with d{6) = 1; 

then T] t is ergodic and its unique equilibrium measure is determined for any y G Y 
by 

u(y)= c (^o=,{^ = ^}' 

d(6) = 1 

where r is the hitting time of {9 G : d{6) = 1} U {p} for £ t . 
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Proof, (a) Since < d(y) < 1 for all y G Y, we can define a Markov chain £ t with state 
space Y = YUjp} and transition probabilities given by 

f eZ(£)P 6=§ {6 = y'} , i/ ii'eY, 



P& =i? {li = y'} := { i-d(y) 



i/ y G Y, £/' = p 



i , if y = y' = p 

Note that the new state p is an absorbing state of the process £ t . 
Define H : X x Y -> R by 



#(x,y) := 




#0,y) , */ y e Y, 



«/ y = P 



(12) 



It is trivial that (2) holds for x G X and p since H(x, p) = 0. On the other hand, 
if y G Y, then 



E 



»?o=a; 



H(r)i,y) = E Vo=x [H(r)i, 



d(y)Eto=y [H{x,^)] 

d(y) H{x, y')Pt =v {6 = y'} 



i = y' 



s/'eYu{p} 



E 



£0=3/ 



Therefore, 774 and £t are dual to one another with respect to H . 
(b) Note first that by i, the set of linear combinations of H(-,y), with y G Y, is 
dense on C(X). 

Denote by := 6 9 : = 1} U {p} the set of all absorbing states of £ f and 
denote by r the hitting time of for £ f , that is, 

r := min < s > : 



js > : i G §} . 



From hypothesis ii, we obtain that there exists a < 1 such that for any non- 
absorbing state y of & we have that d(y) < a and P| 0= ^ ^ £1 ^ @ f < a. It follows that 
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< « E p &h* {£-1 = 2 } 



zeY\e 



Therefore, by induction and noting that P^ )= - {r = oo} < P^ o= - |£ s G" 0| for 
any s > 1, we get Pg 0= £ {t = oo} < a s and, since it holds for all s > 1, we deduce 
P &=s {r < oo} = 1. 

Defining, for any probability measure /i G "P(X) and any y G Y, : = 

J x y)d/i(x), it follows that for any ?/ G Y 



lim/} s (£) = lim E* - 



6»ee 



+ lim E 



T>s P^{r> S } 



6»GG, d(6»)=l 

where the first equality follows from the fact that rjt and £t are dual to one another 
with respect to i7 and equation (10), and the last one is due to in, H(-,p) = and 
= {6 G O : d{9) = l}U{p}. 

□ 

Remark 2.2. From the proof of Theorem 2.1, it is straightforward to check that 
condition ii can be replaced by ^^ 0=y {V < oo} = 1 for all y G Y, where r' t/ie hitting 
time of /or since implies P| 0= ^ {t < oo} = 1 /or all y G Y . 



Remark 2.3. As equation (7) can be seen as a discrete-time counterpart of (4), it is 
interesting to compare the form of the limit measure (11) with the expression given in 
(2.4) of [7] for continuous time: 

%) = E^ =y 



<oo} 
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Note that, in [7], the expression is given through the dual process £ tj studying its 
complete evolution from zero time to the time of absorption, while here it is given only 
through the absorbing point of the modified process £t instead of the original dual process 



Remark 2.4. Theorem 2.1 gives conditions for the ergodicity of the process and, in that 
case, the form of the limit measure. Note that, even when the process is not ergodic, 
equation (7) and the computations made in the proof of the theorem can be used to study 
the behaviour of r\ t as a function of the evolution of £ t . 



3. Nearest-neighbour probabilistic cellular automata 

In this section we will use the duality equation (7) to study conditions for the ergodicity 
and the limit behaviour of some PCA models evolving on Z. 
First, let us introduce a definition: 

Definition 3.1. We say that a PCA rj t with state space X = W z = {1, . . . , M} z is a 
nearest-neighbour PCA if N = {—1, 0, 1}. 

Thus, the evolution of such a process is given for any i,j,k,m G W, z G Z and 
seNby 

F{t] s+1 (z) = m\ 7] s (z - 1) = 2,77,0) =3,Vs{z + 1) = A;} = Pi jk , m - 

We will treat two wide classes of nearest-neighbour PCA, including classical voter 
models, linear and non-linear voter models, noisy voter models, biased voter models and 
some competition models. 



3.1. Multi- opinion voter models 

Throughout this subsection, rjt will be a nearest-neighbour PCA with M states, and the 
state space of the dual processes will be Y = Y M ~ l (where Y is the collection of finite 
subsets of Z). 

Let H : X x Y — > M be the function defined for all x = {x(z) e W} z£ i 6 X and 
A = (A)i<i<A/-i e Y by 

f 1 , i /x W = i ,V,eX„V i = l,...,M-l 
, otherwise 

Our aim is to construct a Markov chain At with state space Y and a function 
d : Y — > [0, oo) that verify equation (7) with H defined in (13). 

Suppose rjt and A t = (Ai it , ■ ■ ■ , Am-i,*) are dual to one another with respect 
to (H,d), for some function d, and let x G X and A = (Aj)i<j<Af-i G Y, where 
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{A;}i< <m_i is a pairwise disjoint family of subsets of Z. Then, equation (7) can be 
written as: 



=x {ni{ z ) = e A u l < i < M - 1} = d(A)P Ao =A {x(z) = i,Vz G A iA ,l < i < M - 1} 

= d(A)P Ao=A c 1 < i < M - 1} , 



(14) 



where R l x := {z G Z : = i}. Note that given x and A, then the left hand side can 
be written as the product 

M-1 

W vo=x { Vl {z) = i, Vz e A, 1 < i < M - 1} = [ F V0=X { Vl (z) = i, Vz e 



M-1 

= n n^o-^i^^o, 

i=l 2gAi 

as long as Al, . . . , Am-i are pairwise disjoint. This factorization, due to the conditional 
independence of different particles, leads us to consider dual chains where the sets 
Ai iS , . . . , A M _i jS evolve independently as long as they are disjoint; that is, given pairwise 
disjoint Ai, . . . , Am-i, we require 

M-1 

Pa () =a {Ai tl = B t : 1<1<M-1}=Y[ P Ao=A {Ai,! = B t } 

i=l 

for all pairwise disjoint B\, . . . , Bm-i G Y. Indeed, we will consider processes in which 
the evolution of points belonging only to one component is independent of the rest. 
More precisely, for pairwise disjoint Ai iS , . . . , Am-i, s , a singleton {z} contained in A m , s is 
substituted by B m (z), a subset of {z — 1, z, z + 1}, independently of the rest of singletons 
in A m>s , with probabilities n m given by 



B m (z) = {z — l,z,z + l} 



B m (z) = {z,z + 1} 




B m (z) = {z-1} 



(15) 



B m (z) = {z-l,z + l} 



{*} 



B m {z) = {z-l,z} B m (z) = {z + 1} 



The set A m , s+ i is the union of the resulting sets from the updating of all {z} C A m , s , 
that is A myS+ ! := [J z€A B m (z). If the sets Ai are not pairwise disjoint, the point 
A = (Ai, . . . , Am-i) is absorbing for the dual process. Note that the evolution of this 
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process can be seen as a multi-type coalescing branching process. Moreover, we will 
take the function d with the form 

M-l 

d(A) = J] d\ A \ (16) 

i=l 

with di G [0, +oo). 

Now, let x G X and AG Y. If A\, . . . , Am-i are not pairwise disjoint, then both 
sides of (14) are zero. If A±, . . . , A M _i are pairwise disjoint, we have, due to the evolution 
of the dual, and noting that Rl, . . . , R^f -1 are pairwise disjoint by definition, 

M-l 

d(A)P Ao=A {A iA C Ri : 1 < i < M - 1} = J] d[ M F Ao=A {A iA C 4} 

i=l 
M-l 

= n^ Ai1 n ^w«)U}. 

i=l {^}CAi 

where Bj(z) is the one-step evolution of {z} C Aj in the dual process. 
Therefore, equation (14) takes the form 

M-l M-l 

n n = i} = n n ^ p a =a qk}- 

i=l zeA t i=l {z}CAi 

Now, we look for conditions on the probabilities ir m such that rjt and At are dual 
with respect to (H, d). We must have f m = x {vi( z ) = m } = ^mIPAo=A {B m (z) C R™} for 
all x G X, A = (Ax, . . . , i4j\f_i) pairwise disjoint, z£Z and m = 1, . . . , M — 1. These 
equations take the following form 



Pijk,m 


~ d 7T 

u "m " m 






Pmjk,m 


= d m [nl + *4J 






Pimk,m 


= d m [irl + 7i c m ] 






Pijm,m 


= d m [nl + ir r m ] 






Pmmk,m 


— A L-0 + TT 1 4- 7T C 


+^} 




Pmjm,m 


= d m [tT^ + 7T^ + -K r m 


+<] 




Pimm,m 


— A \>jr% _|_ «r C _|_ ir r 


+ c] 




Pmmm,m 


— A \Jb + it 1 + 7T C 


1 m 1 m 1 "m 


_|_ n cr _|_ jj-icrl 
1 m 1 m J 



(17) 



with i, j, k E W \ {m}, together with 

7T m > 0, 7T m + 7T^ + 7T C m + 7T^ + 7T,^ + 7T^ + 7T^ + TT^ = 1. (18) 

Note that the last line of (17) is actually 

Pmmm,rn d m . 

The existence of solutions for equations (17) and (18) requires extra assumptions on 
the process r] t . In fact, such equations may only have solutions if for any m G W\ {M} 
each of the parameters Pijk, m ,p m jk, m ,Pimk,m, ■ ■ ■ ,Pimm,m does not depend on k ^ m. 
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That is to say that the probability of a site assuming a state m at time s + 1 is a function 
of the positions of the state m in its neighbourhood at time s. Moreover, when (17) has 
a solution, we can note that p m mm,m = implies Pijk, m = for all k e W, so the 
state m could not be taken by r\ t . Hence, it is reasonable to assume p mm m,m > for all 
meW\ {M}. Thus, < d m < 1 for all m = 1, . . . , M - 1. 

Therefore, for any m G W \ {M} and i, j, fc/mwe can denote p m := Pijk,mi which 
can be interpreted as the probability that r)i(z) spontaneously assumes the state m (i.e. 
when m is not present in the neighbourhood of the point). Thus, p m mm,m is equal to 
the probability that 771(2:) does not spontaneously assume any state k 7^ m, and for 
m — 1, . . . , M — 1, we have 

dm 1 Pmmm,M P ~\~ Pm 

where p := J2n=i Pn- Hence, (16) becomes 

M-l 

d(A) = II (1 - Plll ,M -p + Pi) lM . (19) 

1=1 

The equations (17) and (18) have solutions given by: 
7T - — V 

~i {Pmkk,m Pm) 
" ~j {Pkmk,m Pm) 
= ~1 \Pkkm,m Pm) 

Y (20) 

^"m — ~j (Pmmfc,m H~ ^Pm Pmkk,m Pkmk,m) 

~j (jPrnkm^m T^m — Pmkk.rn Pkkrn,m) 

— : ~3 (jPkmm,m H~ Pkmk.rn Pkkm,m) 

T^rn j \Pmmm.m ~t~ Pmkk.rn i Pkmk,m T" Pkkm,m Pkmm,m Pmkm.m Pmmk.m Pm) 

dm 

with k m, provided that 

Pmkk,m _ Pm 
Pkmk,m _ Pm 
Pkkm,m _ Pm 

Pmmk,m — Pmkk,m Pkmk,m Pm ("^-0 

Pmkm,m _ Pmkk,m Pkkm,m Pm 

Pkmm,m — Pkmk,m Pkkm,m Pm 

Pmmm,m Pmkk,m Pkmk,m Pkkm,m _ Pkmm,m Pmkm,m Pmmk,m Pm 

for all m — 1, . . . , M — 1, k G W \ {m}. The inequalities (21) impose lower and upper 
bounds on the velocity for the growth of the probability that a site assumes the opinion m 
in the time s + 1 because of the occurrence of opinions m in its neighbourhood in time s. 



Thus, we have proved: 
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Theorem 3.2. Let rjt be a nearest-neighbour PC A with state space X, whose transition 
probabilities are such that for any m = 1, . . . , M — 1 each of the parameters 
Pijk,m,Pmjk,m,Pimk,m, ■ ■ ■ ,Pimm, m does not depend on k ^ m and satisfy inequalities 
(21). Then there exist a Markov chain A t with state space Y, H : X x Y — > R, and 
d : Y — > (0, 1], such that rj t and A t are dual to one another with respect to (H,d). 
Furthermore, if p mmm , m = 1» f or all m = 1, . . . , M — 1, then r] t and A t are dual to one 
another with respect to H. 

□ 

The following result follows from the application of Theorem 2.1 and Theorem 3.2. 

Theorem 3.3. Let r\t be a nearest-neighbour PC A whose transition probabilities are such 
that for any m = 1, . . . , M - 1 each of the parameters Pijk,m,Pmjk,m,Pimk,m, ■ ■ ■ ,Pimm,m 
does not depend on i,j, k ^ m and satisfy inequalities (21), and let A t = (A,t)i<i<M-i 
be the process on Y defined by the transition probabilities (15) and (20). Then, the 
process is ergodic if and only if any of the following conditions holds: 

(i) Pm,M > for all i G W \ {M}, that is r] t admits spontaneous changes from any 
pure states i ^ M to the state M; 

(ii) 3m, n G W \ {M}, m ^ n, with p m ,Pn > 0, that is t] t admits spontaneous changes 
to at least two distinct states; 

(Hi) 3m G W \ {M}, with p m > and p mmm ,M > 0; 

(iv) 3m G W \ {M}, with p m > 0, and such that the process A t starting at A = 
(0, . . . , A m , . . . , 0) has probability 1 of extinction for all A m G Y. 

Furthermore, when any of the above conditions holds, the unique equilibrium 
measure for rj t is defined for any A = (Ax, . . . , Am-i) G Y, by 

z>(A)=P Ao=A {A T = (0,...,0)}, 

where A t is the process defined on Y U {p} with transition probabilities given in (12), 
and r is the hitting time of {(0, . . . , 0)} U {p} for A t . 

Proof. We first check the sufficiency of the conditions for ergodicity. 

With the previous construction, r\ t and A t are dual to one another with respect to 
(H, d). Also, by Lemma 3.1 of [7], the set of linear combinations of functions if (•, A) is 
dense in C(X). 

The set of absorbing states of A t contains (0, . . . , 0) together with all A = 
(A±, . . . , Am-i) which are not pairwise disjoint. Moreover, let us denote A the set 
of values m £ { 1 , . . . , M — 1}, if any, such that 

Pm Pmkk,m Pkkm,m Pmkm,m 0, Pkmk,m Pmmk,m Pkmm,m Pmmm,m- 
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Then, if m G A, we have 7r^ = 1 so the points in A m do not change in time. Thus, 
the set of absorbing points of A t also includes the points A = (A±, . . . , Am-i) such 
that At = for all % & A. 

Let us see now that under any of the conditions i, ii, Hi or iv, it holds {8 G O : 
d{6) = 1} = {(0,...,0)}. Obviously, tZ((0, . . . , 0)) = 1. Now, for . . . , not 
pairwise disjoint with Aj fl A,- ^ 0, we have, under any of i, ii, Hi or iv, that at least one 
of di — 1 — Pi^A/ — P + Pi or dj — 1 — Pjjj,M — P + Pj must be smaller than 1, so from 
(19) we get d(A) < 1. Last, let j 6 i; since pi = we have that dj = 1 — pm t M — P', 
note that under i we have Pm t M > 0, and under any of ii, Hi or iv, there must be 
some m E W \ {{M} U A) with p m > 0; therefore under any of i, ii, Hi or iv, we get 
di = 1 — pui ; M — P < 1 — Pm,M — Pm < 1 an d d(A) < 1. Note that in particular, 
{06 0: = 1} = {(0, . . . ,0)} implies that the condition Hi of Theorem 2.1 holds 
with c(0) = 1. 

Suppose z holds, therefore it follows d(A) < max {1 — Pkkk m} < 1 for any 

kew\{M} 

non-empty set A G Y. If ii holds, then d(A) < max{l — pk} < 1 for any non- 

k=m,n 

empty set A G Y. If Hi holds, then for any non-empty set A G Y we have 
d(A) < max{l — p m , 1 — Pmmm,M} < 1- Hence, in all the cases, we apply Theorem 2.1 
to conclude. 

Now, suppose iv holds. We can consider that for any k G W \ {M}, k ^ m, we 
have pk = 0, and also p mm m,M = 0, since if the contrary holds, we could conclude the 
ergodicity by ii or Hi. Hence, for any k ^ m, it follows that irf. = 0, while > 
(see (20)). Notice that Y = Y : U Y 2 , where Y 1 := {B G Y : 3k ^ m, 5 fc ^ 0} and 
Y 2 := {B G Y : VA; ^ m, 5 fc = 0}. For all A G Yi we have d(A) < 1 - p m < 1, and 
if A = A, then A s G Y x for all s < r. Thus, the process A t starting in A G Yj 
has probability 1 of being absorbed in p. On the other hand, for all A G Y 2 we have 
d(A) = 1 and if A t starts in A G Y 2 , then we have A s G Y 2 for all s < r, and it 
coincides with the process At starting in A G Y 2 which has probability 1 of being 
absorbed in (0, . . . , 0). Therefore, in both cases, Pa =a i T < °°} = 1- Thus, from the 
proof of Theorem 2.1 we conclude the ergodicity of r\ t (see Remark 2.2). 

Furthermore, in all the cases we get the expression of the equilibrium measure for 
r) t from Theorem 2.1. 

For necessity we only have to show that the process is not ergodic in the two 
following cases: (a) pi = for alii = 1, . . . , M — 1 and there exists j ^ M such that 
Pjjj,M — 0; or (b) 3!m G {1, . . . , M — 1} with p m > 0, p mmm ,M = 0, and there is some 
A m G Y such that the process At starting in A = (0, . . . , A m , . . . , 0) has probability 
strictly less than 1 of extinction. Case (a) is immediate since the probability measures 
5m and 5j, concentrated on pure state configurations x' G X and x" G X, x'(z) = M 
and x"(z) = j for all z G Z, respectively, are both equilibrium measures for the process. 

For case (b), note first that the measure 8 m is an equilibrium measure for the 
process. Now let x G X with x(z) ^ m for all z G Z and A = (0, . . . , A m , . . . , 0) 
such that At has a probability less than 1 of extinction. Note that the set Y(m) : = 
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{B G Y : Bi = 0, Vi 7^ m} is closed for the evolution of the dual process (that is, 
A t G Y(m) for all t > 0), and d(B) = 1 for all B G Y(m). Then, for the process A t 
enclosed in Y(m), the duality relation (7) takes the form of (2), and for the particular 
value of x and A we have 

IV*. {rjs(z) = m, Vze A rn } = P A()=A {A m , s C R™} 

= Pa„=a {A m , s = 0} 

< F A()=A {r < oo} < 1. 

Then, limsup^^ f Vo=x {Vs(z) — m, Wz G A m } < 1 so the process starting with a 
measure with support in {x G X : 7^ m, Wz G Z} is not converging to <5 m and the 

process is not ergodic. 

□ 

Remark 3.4. Theorems 3.2 and 3.3 can be extended to a wide class of PC A 
with finite range interaction on the integer lattice Z d (1), which are such that 
P {rj s+ i{z) = w\ T) s (z + n) — w n , Vn G N} does not depend on w n 7^ w. In fact, for 
any d and N C Z d we can define a dual process (j on Y = y M ~ 1 j where Y is the 
collection of finite subsets of Z d , and a point z belonging to A m , s is substituted by 
a subset of {z + n : n G N}. Thus, in an analogous way to that in (21) it is possi- 
ble to find inequalities that delimit the class of PC A for which we can extend our results. 

The PCA model presented in this subsection recover many one-dimensional models 
that have been referred in the literature. In fact, our model recovers M-state linear and 
non- linear PCA (which were studied in the 2-state version by [13] and [14]), presenting 
new results for them. 

Example 3.5. (Multi-opinion noisy voter model) Suppose r]t is defined as follows: 
For each m G W = {1,...,M} ; take q m G [0,1], and define q = J2 meW q m ; set 
a, A 7 ^ 0; suc h that a + /9 + 7 = l — q; finally suppose r\ t has transition probabilities 
given by 

Pijk,m ■= al{i}(m) + pl {j} (m) + ^l {k] (m) + q m . 

Note that for any m G W and i,j,k G W \ {m} we have q m = Pijk, m ='■ Pm 
(which here is also defined for m = M). Using the language of voter models, q m is the 
probability of a voter assuming the opinion m, independently of her own opinion, as well 
as the opinions of her neighbours in the previous step, while a, [3 and 7 are the weights 
given by each voter to her own opinion and her neighbours ' opinions. 

Corollary 3.6. Let rj t be a multi-opinion noisy voter model. Then there exists a 
coalescing random walk A t on the family of all finite subsets of Z such that rjt and 
A t are dual to one another with respect to (H,d) defined in (13) and (19). Moreover, 
rjt is ergodic if and only if there exists m G {1, . . . , M} such that p m > 0. 
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Proof. It is easy to check that for any m = 1, . . . , M — 1 the parameters Pijk, m ,Pmjk,mi 
Pimk,m, ■ ■ ■ ,Pimm,m do not depend on i, j, k ^ m and satisfy inequalities (21). Hence, we 
can apply Theorem 3.2 to deduce the existence of a dual with respect to (H, d), where 
H is given by (13) and d is given by (19). 

For the second part of the result, if p m = for all m, then any probability measure 
concentrated on a pure state configuration x G X, x(z) = m for all m G Z and any 
m G W, is an equilibrium measure. For sufficiency, note that (21) holds regardless of 
the order of the opinions, and thus we can set M to have Pm > and ergodicity follows 
from Theorem 3.3 i. 

□ 

Example 3.7. (Domany-Kinzel model) If we consider the case W = {0, 1} and 
such that the probability of a site z assuming 1 at time s + 1 is a function of the number 
of Is in the neighbourhood {z — 1, z + 1} at time s, then we have the Domany-Kinzel 
model (6) with parameters 

a o '■= Pooo.i = Poio,ij 

a i '■= Pioo.i = Piio,i = Pon,i — Pooi,ii 

a-2 ■= Pm,i =Pni,i- 

We recall that the Domany-Kinzel model was introduced in [15] and has been ex- 
tensively studied (see for instance [12], [16], [17] and [18]) due to its useful applicability 
in percolation theory and phase-transition theory among others. 

Corollary 3.8. Let rjt be the Domany-Kinzel model (6) with parameters < oq < ai < 
a 2 < 1 . Then there exists a coalescing random walk A t on the family of all finite subsets 
ofTL, such that rj t and A t are dual to one another with respect to the pair (H, d) defined 
by (13) and (19). Furthermore, if one of the following conditions holds 

(i) «o + a 2 > 2ai and a-i < 1; 

(ii) ao + a2 < 2ai and a > 0; 
(Hi) oq > and < 1; 

(iv) a = 0, ai < 1/2 and < I; 

(v) ao > 0, ai > 1/2 and a 2 — 1; 

then rjt is ergodic. 

Proof. If ai < o( a o + a 2), then r\t satisfies all hypotheses of Theorem 3.2 (identifying 
the state here with the the state M in Theorem 3.2). Hence, we have that A t is 
defined by the transition probabilities 7r := ao/a2, ir e = ir r := (ai — ao)/a2, and 
ir er := (a 2 - 2a x + a )/a 2 - 
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On the other hand, if a% > |(d + d 2 ), then we can consider the description of the 
process by the probabilities of the occurrence of 0, that is the Domany-Kinzel model 
with parameters 

a* := 1 - a 2 _i, i — 0, 1, 2. 

It is easy to see that < < a* < a\ < 1, and d^ < ^(dg + d 2 ). Therefore, 
we can again apply Theorem 3.2 (now the state 1 plays the role of the state M in 
Theorem 3.2) and find A^, which is defined by the transition probabilities 7rf := d^/d^, 
*i = < := (d^ - d^/d^, and vrf := (a* - 2aJ + a* )/a* 2 . 

Furthermore, we can obtain sufficient conditions for ergodicity: 

(i) Identifying the state M with in Theorem 3.2 and pm ; o = 1 — d 2 > 0, the result 
follows from Theorem 3.3 i. 

(ii) We now identify M with 1 so pooo,i = 1 — d 2 = do > and the conclusion follows 
from Theorem 3.3 i. 

(iii) We identify M with or 1 depending on (a 2 + ao)/2 being greater or smaller than 
a±. In both cases, for i^Mwe have pm t M > and the result follows from Theorem 
3.3 i. 

(iv) Suppose ao = 0, di < 1/2, and a 2 < 1. We have two possible cases: 

(a) If di < d 2 /2, then we can use the dual process A t and, since pm,o = l — a 2 > 
and the state plays the role of the state M in the Theorem 3.2, we can deduce 
the ergodicity of r]t from Theorem 3.3. i. 

(b) If di > d 2 /2, then consider the dual process AJ 5 (that is, with M — 1, pooo.i — 0, 
Po = 1 — a 2 > 0), starting from any Ao G Y. Let us show that A£ has 
probability 1 of extinction. 

Consider the branching process Xt on Z defined as follows: Let \s be the 
number of individuals in the sth generation of some population; assume each 
individual independently will give rise to offsprings with probability 7rf, 1 
offspring with probability (7rf + 7r£), or 2 offsprings with probability 7rf"; then 
Xs+i is the number of individuals in the next generation, that is the total 
number of offsprings generated by the individuals in the generation s. Since 
di < 1/2, it follows that the expected number of offsprings of each individual 
is (nl + 7r£) + 2ir l J = 2di < 1, which implies the branching process Xt vanishes 
with probability 1 (see [19]). Identifying xt starting with |A | individuals as 
an upper bound on the number of individuals of A^, we have that the dual 
process A^ has probability 1 of being absorbed in and we can use Theorem 
3.3. iv to deduce the ergodicity of rj t . 

(v) If d > 0, di > 1/2 and d 2 = 1, then Oq = 0, a\ < 1/2 and a 2 < 1, and therefore we 
get the result reasoning as in iv. 



□ 
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Remark 3.9. As we said in the Introduction, the inclusion of the function d in our 
definition of duality (7) allows, in some cases, to analyze processes which do not have a 
dual with respect to a given function H in terms of the classical duality equation (2). In 
fact, from the construction of Theorem 3.2, we can see that if we consider the classical 
duality (2) with respect to the function H defined in (13) we would obtain a system of 
equations similar to (17), but without the factor d m multiplying the right side of the 
equations. Therefore, for the existence of a solution of the system, we would need the 
extra hypothesis p m mm,m = 1 for all m = 1, . . . , M — 1. In particular, M > 3 would 
imply that every state m = 1, . . . , M — 1 should have p m = and p mm m,M = 0; while if 
M = 2, we could have p\ > but not pin,2 > 0. Note that the additional hypothesis of 
Pmmm,m — 1 f or all m ^ M would imply that if M > 3, then rj t is not ergodic; and if 
M = 2, then r] t would be ergodic if and only if p\ was large enough in order for condition 
of Theorem 3. 3 iii to hold. 

In Example 3.5, if we had used the classical duality equation (2), we would be 
restricted to consider models with every q m = (if M > 3) or with only one q m ^ (if 
M = 2). On the other hand, in Example 3.7, and due to symmetry of states '0' and l l ' 
in the Domany-Kinzel model, we could have either ao > or 02 > 0, but not both ao 
and 02 positive. 



3.2. Multi-state monotone models 

In the previous subsection we showed that a necessary condition for a process to have 
a dual with respect to (H, d) defined in (13) and (19) is that for any m — 1, . . . , M — 1 
each of the parameters p..., m depends only on the presence (and position) of the value 
m in the neighbourhood of the point. However, there are models (for instance, when 
the set W has an order) where not only the presence (and position) of the value m is 
important for the probability of changing to m but also the presence (and position) of 
values bigger or smaller than m is. In order to find a dual for these processes, we will 
use a duality function H different to (13) and allow more transitions for the dual than 
we did in Subsection 3.1; as we will see, in this case, the evolution of any two subsets 
Ai and Aj has a greater interdependence than in the non-biased case. 

For simplicity, we consider only the case of 7] t being a nearest-neighbour PCA for 
which the probability of a fixed site assuming some state at time s + 1 does not depend 
on its own state at time s. That is, we assume that 

Pik,ra ■= Pijk,m, Vi, j, k,mEW. (22) 

Throughout this subsection, the state space of the dual process will be Y = 
{A G Y M ~ 1 : Ai C ■ ■ ■ C A M -i) and the function i/:XxY^I, defined by 

(l , ifx{z)<k, VzeA k , Vk=l,...,M-l 
, otherwise 

The class of dual processes we consider is composed of processes A t (= £ t ) starting 
at A G Y whose transition probabilities from A s = (A itS )i<i<M-i to A s+ i = 
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(Ai tS+ i)x<i<M-i are defined as follows. Let (0, . . . , 0) be an absorbing state for A t and, 
under the convention that A ^ s = 0, assume that for any k,m,n £ {1, . . . , M — 1} each 
singleton {z} £ Ak, s \Ak-i tS will generate an element independently of other z' £ Am-i, s , 
B(z) = (Bi(z))i<i<M-i £ Y as follows: 

• with probability nf, 

Bi(z) = Q, \<i<M-l; 

• with probability 7r| m 



Bi(z) 

with probability n r k m 



BAz) 



, if 1 <i <m 

{z-1} , if m < i < M - 1; 

, if 1 < i < m 

{z + 1} , if m <i < M- 1; 



with probability 7r fcj , 



{,-1} 
{*+!} 



if 1 <i < min{m,n} 

if m < i < n 
if n < i < m 



{z — 1, z + 1} , max{m, n} < i < M — 1; 



with 



M-l 



M-l 



M-l M-l 



n k + Yl n l™ + Y n k,n + Y Y nk > mn 



m=l n=l m=l n=l 

We define A s+1 = (A, s +i)i<;<m-i, where 

A i>s+1 = y ^(2). 

zeAif-1 

Note that A s+1 also belongs to Y. 

Equation (7) is, in this case, for each x £ X, A £ Y, 

F m=x { m (z) < k : Vz e A fcl 1 < k < M - 1} = d(A)P Ao=A {A fe .i Cft*:l<fc<M-l} 

where 72.* := {2 £ Z : x(z) < k}. 

The left hand side of (24) takes the form: 

M-l 

n n v vo=x { m (z)<k} 

k=l zeA k \A k ^ 



(24) 
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M-l 



d(A) =][d 



A fe \A-il 



fc=l 



for some di, . . . , du-i G [0, 1] (the product over the empty set is taken as 1), the right 
hand side of (24) is 



M-l 

II II 4P Ao =A{£ m (z)C7C:l<m<M-l}. 

k=l {z}QA k \A k _ x 

Then, we must find solutions for the set of equations 

P w =. 07i (*) <k} = d k ¥ Ao=A {B m (z) C : 1 < m < M - 1} , (25) 

for all £ G X, A G Y, fc = 1, . . . , M - 1 and 2 G A k \ A k - X . Writing S* = ElUi P*i,m 
and noting that, for 1 < z, j < M, x(z — 1) = i, x(z + 1) = j, we have 

f 

[ {z-l,z + l} 
The equations (25) are 



if i < m < j 
if j < m < i . 
for m> i,j 



qk 

J M,M 



ok 

J M-1,M 
Qk 

°M,M-1 
ok 

J M-l,M-l 



ok 

°M-2,M-1 



dk4 



d k K + ^Im-i] 
dk [4 + 7Tj |W _J 

d k K + ^k.M-l + ^.M-l + TTfc.Af-1 M-l] 

4 [vrfj + 7T^ M _ 1 + 7rf )M _ 2 + VT^ M „ 1 + 7T fcj M-l M-l + 7Tfc,M-2 M-l] 



(26) 



M-l 



7l 



M-l 
n=j 



M-l M-l 

EE 

m=i n=j 



#1 



M-l 



M-l 



M-l M-l 



»i = l 



71=1 



m=l n=l 
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for k = 1, . . . , M — 1. Due to the last line of the above equations, we deduce that 
dk = S^i and then 

M-l / k \ l A kW-ll 

fc=l \m=l / 

We assume pn t i = 5^ > 0, for otherwise, under (26), Sjj = for all i, j G {1, . . . , M} 
and opinion 1 cannot be taken by the process. Therefore, dk = > for all 
k G {1, . . . , M — 1} and the solution of (26) is given by: 



7T 



n k,m 



k,m 



qk 

°M,M 

dk 



ok _ qk 
°m,M °m+l,M 

d k 



qk _ qk 

°M,m °M,m+l 



d 



k 



qk qk qk _i_ qk 



dk 



since 



M-l M-l 

} M,M 



dk ^k,m — {Sm,M ^m+l,Af) — &i,M 



M-l M-l 

ffc 

'M,M> 



j \ ^ r \ ^ / qk _ qk \ qk _ qk 

k 2.^1 k > n ~ 2.^1 \ M,n °M,n+l) — J M,j °A 



M-l M-l M-l M-l 



dk ^k,mn ~ {^m,n ^m+l,n ^m,n+l + ^m+l,n+l) 



m=i n=j m=i n=j 

M-l 



(2* 



(29) 



Ef qk qk _ qk _j_ qk \ 

i^mj ~~ °m,M ~ °m+l,j ' °m+l,M) 

m=i 

cfc cfc qk I cfc 

_ D i,j °M,j °i,M °M,M- 

In order to have the nonnegativity of (28), we need 

S m ,M — Sm+l,Mi ^M,n — ^M,n+li (30) 

Sm,n ~ ^m+l,n ~ ^m,n+l + ^m+l,n+l — 0' (31) 
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for all 1 < m,n < M - 1. Note that (30) and (31) together imply > S^ n for all 
1 <i <m < M, I < j <n < M and k = 1, . . . , M, that is 

k k 

Pij j > ^2 Pmn ' e 

which means that the process r\t is monotone with respect to the coordinatewise order 
in X induced by the total order l<2<---<Minjy (see [21]). Therefore, condition 
(30) can be substituted by the monotonicity of the process. 

On the other hand, conditions (30) and (31) can be written equivalently as saying 
that Sfj = 1 — F k (i,j) for each k = 1, . . . M — 1, where F k is a distribution function on 
R 2 . 

We have proved the following result: 

Theorem 3.10. Let rjt be a nearest-neighbour PCA with state space X whose transition 
probabilities satisfy (22). Then, there exist a dual process A t with state space Y such 
that rjt and A t are dual to one another with respect to (H, d) defined in (23) and (27) if 
and only if rjt is monotone and its transition probabilities satisfy (31). Furthermore, if 
Pn i = 1, then r] t and A t are dual to one another with respect to H . 

□ 

Corollary 3.11. A monotone nearest-neighbour PCA whose transition probabilities 
satisfy (31) is ergodic if Pn,M > 0. In such case, the unique equilibrium measure for r\ t 
is defined for any A = (A\, . . . , Am-i) € Y } by 

^(A)=P Ao= a{a. = (0,...,0)}, 

where A t is the process defined on Y U {p} with transition probabilities given in (28), 
and t is the hitting time of {(0, . . . , 0)} U {p} for A t . 

Proof. By Theorem 3.10 there exist a Markov chain A t and a pair of functions (H,d), 
such that rjt and A t are dual to one another with respect to (H, d). 

We now apply Theorem 2.1 to get the result. In fact, by Lemma 3.2 of [7] the set 
of linear combinations of the functions H(-,A) with A e Y is dense in C(X). Also, 
(0, . . . , 0) is the unique absorbing state 9 of A t such that d{9) = 1, and H(-, (0, . . . , 0)) = 
1. Moreover, for any A ^ (0, . . . , 0), d(A) = Ul=i\^-EjLk + iPii,j) lAkXAk - l{ < 1~Pii,m- 
Thus, if pn t m > 0, hypotheses ii and Hi (with c(0) = 1) of Theorem 2.1 hold and the 
result is proved. 

□ 

Remark 3.12. Note that as in the case of the multi-opinion voter models, we can 
extend Theorem 3.10 and Corollary 3.11 to a wide class of multi- state monotone biased 
models with state space W zd and finite range interaction. 
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We notice that the processes r\ t that verify the conditions of Theorem 3.10 include 
many PCA models, such as multi-opinion noisy biased voter models and multi-species 
competition models. 

Example 3.13. (Multi-species competition model) Suppose M distinct species 
live in the space Z. Each site of the space contains only one individual. The 
spread mechanism for all species is the same: each individual produces one offspring 
(independently of other individuals) that will leave to colonize a neighbouring site, killing 
the old particle that occupied the site and competing against the offspring of others 
individuals. Moreover, we include the possibility of a spontaneous appearance of an 
individual of any type (for instance, coming from outside). 

We consider such a model described as follows: Suppose r] t is a nearest-neighbour 
process with M states such that the probability Pij jTn that a site assumes the state m in 
the next step when its left neighbour is in the state i and its right neighbour is in the 
state j is given for any j,k ^ m and i < m by 

Pjk,m Pm 

Pim,m Pm (1 P)^m 

Pim,i =Pi + (l-p)(l- a m ) 

Pmi,m Pm (1 P)Pm 
Pmi,i =Pi+ (1 - Pm) 

Pmm,m 1 P ~t~ Pmi 

where p = Y^ m =iPm, & M > oi M -i ■ ■ ■ > oc 2 > 1/2 and f3 M > Pm-i ■ ■ ■ > P2 > 1/2. 

The parameters a m and j3 m can be interpreted as the competitiveness of species m 
with relation to species i < m, and p m corresponds to probabilities that spontaneous 
mutations to m occur. Note that the order in the as and /3s imply that species are 
ordered from weakest 1 to strongest M . Usually, we have a m = (3 m for all m, which 
means there is not a directional component in the competition. 

Competition models belong to the class of biased voter models. The biased voter 
model was introduced in the 2-state version by Williams and Bjerknes [20] as a model 
of skin cancer, and has been applied to study several phenomena (for instance, see [22], 
[23], and [24]). 

Since 

k 

S tj = ^Pe + i 1 ~ P) (C 1 - Ujftii&Kj} + i 1 ~ A)l{j<fc<i} + !{fc>ij}) , 

it is not hard to check that rj t is monotone and (31) holds for it. Therefore, by Corollary 
3.11, if rjt admits a spontaneous change to the state M, then it is ergodic. 
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Example 3.14. Consider the case of monotone processes for which Sf, = g k {i + j), 
where g k : {2, . . . , 2M} — > [0, 1]. In such case, writing I := i + j , condition (31) takes 
the form 

l(9k(l)+g k (l + 2))>g k (l + l), 

which means that (31) is satisfied if and only if the functions g k are convex. In this 
case, we have that rjt is ergodic if 5 , m-i(2) < 1. In fact, if <7m-i(2) < 1, then 
Pn,M = 1 _ •S'li" 1 = 1 ~~ <?Af-i(2) > and we can use Corollary 3.11 to deduce the 
ergodicity of the process. 

Remark 3.15. As in Subsection 3.1, our definition (7) of duality with respect to (H,d) 
allows us to include processes that do not have a dual in terms of the classical definition 
(2). In fact, we could not obtain a version of Theorem 3.10 for classical duality (2) with 
respect to the function H defined in (23) unless we imposed another restrictive condition 
on r] t since we would obtain a system of equations similar to (26), but without the factor 
dk multiplying the right side of each equation. Therefore, we should have S^ — l for any 
k — 1 . . . , M — 1, which is equivalent to imposing the condition pu \ = 1. In such case, 
we would have pn,M = and it would not be possible to obtain a sufficient condition for 
ergodicity in the form of Corollary 3.11. 

In particular, in Example 3.13, classical duality (2) with respect to H would need 
Pm = for any m ^ 1 (which implies pu = 0, and therefore we could not deduce the 
ergodicity), while in Example 3.14 we should impose the condition <?i(2) = 1 (which 
implies <7m-i(2) =1) and therefore we could not deduce the ergodicity either. 

4. Concluding remarks 

In this work, we have presented a new duality equation (7) and shown its usefulness in 
the study of PCA. We believe that there is still room for future research on this topic. 
In particular we are considering the following issues: 

(i) Studying in more detail the equilibrium measures (proportion of values, spatial 
correlations, etc.) when the processes are ergodic in Examples 3.5, 3.7 and 3.13. 

(ii) In the case of processes having a dual with respect to (H, d) in Theorems 3.2 and 
3.10 but not being ergodic, studying their asymptotic behaviour using their dual 
processes. 

(iii) Studying other classes of functions H giving conditions on the transition 
probabilities for the existence of a dual with respect to (H, d). 

Acknowledgments 

We thank the referees for their careful reading and useful comments. 

This work was supported by projects MTM2004-01175 and MTM2007-63769 of 
MEC, and research group "Modelos Estocasticos" (DGA). M. Sobottka was partially 



Dualities for Probabilistic Cellular Automata 



21 



supported by Nucleus Millennium Information and Randomness ICM P04-069-F and 
project DIUC 207.013.030-1.0. 

References 

[1] Lebowitz J. L., Maes C. and Speer E. R., Statistical mechanics of probabilistic cellular automata, 

1990 J. of Statistical Physics 59 117-170. 
[2] Toom A. L., Vasilyev N. B., Stavskaya O. N., Mityushin L. G., Kurdyumov G. L. and Pirogov 

S. A., Discrete local Markov systems, 1990 in:i?. L. Dobrushin, V. I. Kryukov, A. L. Toom, eds., 

Stochastic Cellular Systems: Ergodicity, Memory, Morphogenesis 1-182. 
[3] Liggett T. M., Interacting particle systems, 1985 Springer- Verlag, New York. 
[4] Durrett R., Lectures Notes on Particle Systems and Percolation, 1988 Wadsworth, Belmont, CA. 
[5] Gray L., Duality for general attractive spin systems with applications in one dimension, 1986 Ann. 

Prob. 14 371-396. 

[6] Katori M., Reformulation of Gray's duality for attractive spin systems and its applications, 1994 
J. Phys. A 27 3191-3211. 

[7] Lopez F. J. and Sanz G., Duality for general interacting particle systems, 2000 Markov Processes 

and Relat. Fields 6 305-328. 
[8] Katori M., Konno N., Sudbury A. and Tanemura H., Dualities for the Domany-Kinzel model, 2004 

J. of Theoretical Probability 17 131-144. 
[9] Sudbury A. and Lloyd P., Quantum operators in classical probability theory. II. The concept of 

duality in interacting particle systems, 1995 Ann. Prob. 23 1816-1830. 
[10] Konno N., Dualities for a class of finite range probabilistic cellular automata in one dimension, 

2002 J. of Statistical Physics 106 915-922. 
[11] Konno N. and Ma X., Site-bond representation and self-duality for totalistic probabilistic cellular 

automata in one dimension, 2005 J. of Statistical Mechanics: Theory and Experiment P03002. 
[12] Konno N., Self-duality for multi-state probabilistic cellular automata with finite range interactions, 

2002 J. of Statistical Physics 106 923-930. 
[13] Cox J. T. and Durrett R., Nonlinear voter model, 1991 In: R. Durrett and H. Kesten (eds.), 

Random Walks, "Brownian Motion and Interacting Particle Systems'' Birkhauser 189-201. 
[14] Liggett T. M., Coexistence in threshold voter models, 1994 Ann. Prob. 22 764-802. 
[15] Domany E. and Kinzcl W., Equivalence of cellular automata to Ising models and direct percolation, 

1984 Phys. Rev. Lett. 53 311-314. 
[16] Atman A. P. F. and Morcira J. G., Growth exponent in the Domany-Kinzel cellular automaton, 

2000 Eur. Phys. J. B. 16 501-505. 
[17] Kemper A., Schadschncidcr A. and Zittartz J., Transfer-matrix DMRG for stochastic models: the 

Domany-Kinzel cellular automaton, 2001 J. Phys. A: Math. Gen. 34 279-280. 
[18] Riegcr H., Schadschncidcr A. and Schrcckcnbcrg M., Reentrant behaviour in the Domany-Kinzel 

cellular automaton, 1994 J. Phys. A 27 423-430. 
[19] Athreya K. B. and Ncy P. E., Branching processes, 1972 Springer- Verlag, New York. 
[20] Williams T. and Bjcrkncs R., Stochastic model for abnormal clone spread through epithelial layer, 

1972 Nature 236 19-22. 

[21] Lopez F. J. and Sanz G., Stochastic comparisons for general probabilistic cellular automata, 2000 

Stat. Prob. Lett. 46 401-410. 
[22] Jimenez A., Tiampo K. F. and Posadas A. M., An Ising model for earthquake dynamics, 2007 

Nonlinear Processes in Geophysics 14 5-15. 
[23] Kordzakhia G. and Lalley S. P., A two-species competition model on Z d , 2005 Stochastic Processes 

and their Applications 115 781-796. 
[24] Koulis T. Stochastic modeling of sea ice, 2003 Presented at JSM 2003, San Francisco, USA. 



